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Abstract. We present the lACOB grid-based automatic tool for the quantitative spectroscopic 
analysis of O-stars. The tool consists of an extensive grid of FASTWIND models, and a variety 
of programs implemented in IDL to handle the observations, perform the automatic analysis, 
and visualize the results. The tool provides a fast and objective way to determine the stellar 
parameters and the associated uncertainties of large samples of 0-type stars within a reasonable 
computational time. 



1. Introduction 

The era of large spectroscopic surveys of massive stars has aheady begun, providing us with a 
huge amount of high-quahty spectra of Galactic and extragalactic 0-type stars. The lACOB 
spectroscopic survey of Northern Galactic OB-stars is one of them. This long-term observational 
project is aimed at building a multi-epoch, homogeneous spectroscopic database of high- 
resolution, high signal-to-noise ratio spectra of Galactic bright OB-star^j. Associated with this 
spectroscopic dataset are several working packages aimed at its scientific exploitation. Within 
the framework of the Working Package WPS: Quantitative Spectroscopic Analysis, we have 
developed a powerful tool for the automatic analysis of optical spectra of 0-type stars. The 
tool provides a fast and objective way to determine the stellar parameters and the associated 
uncertainties of large samples of O-stars within a reasonable computational time. Initially 
developed to be used for the analysis of spectra of 0-type stars from the lACOB spectroscopic 
database, the tool is now also being applied in the context of the VLT-FLAMES Tarantula 
Survey (VFTS) project [3], and other studies of stars of this type. 



2. The lACOB grid-based tool 

Apart from the already mentioned characteristics (automatic, objective, and fast), we also aimed 
at a tool that is portable, versatile, adaptable, extensible, and easy to use. As shown throughout 
the following text, this philosophy has guided the whole development of the tool. 

^ Some notes on the present status of lACOB spectroscopic database, the JACOB project, and the various working 
packages can be found in [1] and [2]. 



One of the key-stones of any automatic quantitative spectroscopic analysis is the computa- 
tion of large samples of synthetic spectra (using a stellar atmosphere code), to be compared 
with the observed spectrum. In contrast to other possible alternatives (the genetic algorithm 
-GA- employed by [4], or the principal component analysis -PCA- approach followed by [5]), 
we decided to base our automatic tool on an extensive, precomputed grid of stellar atmosphere 
models and a line-profile fitting technique (i.e., a grid-based approaclo) This option was possi- 
ble thanks to the fast performance of the FASTWIND cod^, and the availability of a cluster 
of ~ 300 computers at the Instituto de Astrofi'sica de Canarias (lAC) connected through the 
condoeEI workload management system. 

The grid is complemented by a variety of programs implemented in IDL, to handle the ob- 
servations, perform the automatic analysis, and visualize the results. The IDL-package has been 
built in a modular way, allowing the user, e.g., to easily modify how the mean values and uncer- 
tainties of the considered parameters are computed, or which evolutionary tracks will be used 
to estimate the evolutionary masses. 

2.1. The lACOB FASTWIND grid 

In the following we outline the main characteristics of the FASTWIND grid used for the analysis 
of Galactic 0-type stars that is presently incorporated within the I AC OB grid-based tool. 

Input stellar parameters considered in the FASTWIND models: 

• Effective temperature (Tgg) and gravity {logg): Table [T] indicates the ranges of T^q 
and log g considered in the grid. Basically, the grid points were selected to cover the region 
of the log Teff - log g diagram where the 0-type stars are located. 

• Helium abundance (ehc) and microturbulence (Ct)- The grid includes six values of 
helium abundance — ehb = Nhc/(Nh+Nhc) — , indicated in Table [TJ For all models, a 
microturbulence Ct,mod = 15kms^^ was adopted in the computation of the atmospheric 
structure, and four values of the microturbulence were considered in the formal solution. 

• Radius {R): Computing a FASTWIND model requires an input value for the radius. This 
radius has to be close to the actual one (which will be derived in the final step of the 
analysis, and hence is not know from the beginning on). We assumed a radius for each 
(Teff, logg) pair following the calibration by [8]. The grid is hence divided in 20 regions in 
which a different radius is considered. For example, for the case of log 5' = 4.0 and 3.5 dex, 
the radii range from 7 to 12 Rq and from 19 to 22 Rq, respectively. 

• Wind parameters (M, Uoo, /?): As conventional for grid computations of stellar 
atmosphere models for the optical analysis of 0-stars, the mass loss rate (M) and terminal 
velocity (uoo) of the wind have been compressed, together with the radius, into the wind 
strength parameter (or optical depth invariant), Q = M /{vooR)^'^ (see [7]). Ten logQ- 
planes were considered for the grid. For each FASTWIND model, M and Voo need to be 
specified. First, a terminal velocity Voo = 2.65fesc was adopted, following [9]; then a mass 
loss rate was computed for the given logQ, Voo, and R. Finally, the exponent of the velocity 
law, /3, was assumed as a free parameter ranging from 0.8 to 1.8. 

^ A similar approach has been recently applied by [6] for the analysis of a sample of 12 low resolution spectra of 
B supergiants in NGC55. 

Among the various available modern stellar atmosphere codes incorporating wind and line- blanketing effects, 
FASTWIND [7] is the fastest. 
* http://www.cs.wisc.edu/condor/ 
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Figure 1. Location of (Teff, logg)- 
pairs considered in the FASTWIND 
grid incorporated into the present 
version of the lACOB grid-based 
tool. Evolutionary tracks from [11]. 
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Table 1. Range of values considered for the six free parameters in the FASTWIND grid, 

logg: [2.6-4.3] dex (with step size 0.1 dex) 

Teff : > 25000 K (step size 500 K), upper hmit defined by the 120 M0 track (see Fig. ^ 

Ct '■ Model: 15 kms~^, Formal solution: 5, 10, 15, 20 kms~^ 

eHe : 0.06, 0.09, 0.13, 0.17, 0.20, 0.23 

logQ : -15.0, -14.0, -13.5, -13.0, -12.7, -12.5, -12.3, -12.1, -11.9, -11.7 

f3: 0.8,1.0,1.2,1.5,1.8 



• Metallicity: A solar metallicity (following [10]) was assumed for the whole gri 

Synthetic lines: The following (optical) lines where synthesized in the formal solution: H^-e, 
HeiAA4026, 4120, 4143, 4387, 4471, 4713, 4922, 5015, 5048, 5875, 6678, and HeiiAA4026, 4200, 
4541, 4686, 5411, 6406, 6527, 6683. 



In order to optimize the size and the read-out time of the grid, only part of the output from the 
FASTWIND models is kept and stored in IDL xdr-files. This includes the input parameters, the 
H/He line profiles and equivalent widths, the information about the stellar atmosphere structure 
and the emergent flux distribution, and the synthetic photometry resulting from the computed 
emergent flux distribution. Using these xdr-files, the size of the grid can be reduced to a 10 % 
of the original one. IDL can restore each xdr-file and compute the corresponding Xj, quantity 
(see Sect. 12. 2p in 0.02-0.1 s per model, i.e., the tool can pass through 80,000 models in 30 
min-1 hour. Following the main philosophy of the lACOB grid-tool, only the reduced grid is 
used by the tool, whilst the original grid is safely kept on hard-disk at the lAC. This allows 
for an easy transfer of the grid to an external disk, hence satisfying our constraint for portability. 

The FASTWIND grid presently incorporated into the lACOB grid-based tool consists of 
~ 192,000 models (~ 32,000 models per He-plane). The reduced size is ~ 34 Gb. This predefined 
grid can be easily updated and/or extended if necessary using appropriate scripts implemented 
in IDL, CONDOR, and LINUX. For example, a new He-plane can be computed and prepared 
for use in ~ 3 days. 



^ We also calculated similar grids for other metalicities {Z = 0.5, 0.4, and 0.2 Zq). 



2.2. The lACOB grid-based tool IDL package 

The guideline of the automatic analysis is based on standard techniques applied within the 
quantitative spectroscopic analysis of 0-stars using optical H/Hc lines that have been described 
elsewhere (e.g. [12], [13]). The whole spectroscopic analysis is performed by means of a variety 
of IDL programs following the steps indicated below. In brief, once the observed spectrum is 
processed, the tool obtains the quantity Xj. (i.e. an estimation of the goodness of fit) for each 
model within a subgrid of models selected from the global grid, and determines the stellar pa- 
rameters and their associated uncertainties by interpreting the obtained Xj, distributions. 

Step 1: The input 

In this first step, the user has to provide the observed spectrum, to indicate the corresponding 

resolving power (i?), and to give some pre-determined information about the star - the projected 
rotational velocity {vslni), the size of the extra line-broadening (Grt), and the absolute visual 
magnitude (My). 

Concerning the grid, the user must select the appropriate metallicity, and indicate the range of 
values for the various free parameters (defining the subgrid of models to be considered within 
the analysis). This latter option allows the tool to be faster (using optimized ranges for the 
various free parameters), or the user to perform a preliminary quick analysis by fixing some of 
the parameters. For example, one can obtain a quick estimate on Teg, log 5, and logQ in less 
than 5 min, by fixing the other three parameters ((^t ^He, and/or j3) . 

Finally, the H/He lines which shall be considered in the analysis and the corresponding weights 
need to be specified. 

Step 2: Processing of the observed spectrum 

In many cases, the observed spectra need to be processed before launching the automatic anal- 
ysis, because of, e.g., nebular contamination affecting the cores of the H and He I lines, the need 
to improve the normalization of the continuum adjacent to the line, and the presence of cosmic 
rays. 

The possibility to correct the observed spectrum for these effects has been incorporated into 
corresponding IDL procedures. The implemented options include (for each of the considered 
lines): local renormalization, selection of the wavelength range of the line, clipping/restoring 
of certain parts of the line, and computation of the signal-to-noise ratio (S/N) in the adjacent 
continuum. 

Once finalized, the processed spectrum is stored. This spectrum and all associated information 
is used in Step 3, and also each time the analysis of the same star/spectrum is (re-)launched. 
This way the user can easily reconsider his/her decisions after a first analysis. 
In this step, the user can select a model from the grid to help him/her with the processing of 
the observed spectrum. One can, for example, make a quick preliminary analysis (see above) 
fixing some parameters and including only few lines, and then use the resulting best model to 
check the processed lines and find constraints for the processing of the others. 

Step 3: Estimation of the goodness of fit for each considered model 



This step, together with Step 4, are the most important ones, constituting the core of the 
program. In the present version of the lACOB grid-based tool we have considered a proce- 
dure as described below; however, we are aware that this procedure can be subject to discus- 
sion/improvements. Having this in mind, and following the philosophy of the tool (regarding 
versatility and adaptability), the corresponding IDL modules have been implemented with the 



possibility to be easily modified. The fast performance of the lACOB grid-based tool makes it 
very powerful to investigate various alternative strategies. 



In our present version, the tool computes, for every model in the subgrid, the quantity 
for each considered line 

""'-N-^h — ^1 — 

where F^ x ™d F^ x are the normalized fluxes corresponding to the synthetic and observed 
spectrum, respectively; aL = (S/N)~^ accounts for the S/N of the line; and Nx is the number of 
frequency points in the line. Under ideal conditions (e.g., for a perfect model, but see below), 
Xj^ corresponds to a reduced x^- 

In a second step, the Xj^ values for each model are corrected for possible deficiencies in the 
synthetic lines (due, for example, to deficiencies in the model, an incorrect characterization 
of the noise of the line, bad placement of the continuum, or bad characterization of the line- 
broadening). To this aim we compute, for each line, the standard deviation of the residuals 
<Tc = o'{Fm,\- Fo^x) from that model that results in the minimum X£ for the given line. Then, 
the following correction is applied: 



2 



Xl, = Xi-^ (2) 

Using these individual X£ ^ values and the weights assumed for each of the considered lines {wl, 
e.g., [4]), a global Xj^ is obtained: 

Nl 

X^ = Y.wlXI, (3) 

L=l 

where A'^^;, is the number of lines. For a large number of frequency points per line, Xj^ should be 
normally distributed. 

Thus, Step 3 provides the values of the quantity Xj^ associated with each of the models in- 
cluded in the considered subgrid. As indicated in Section 12. H this step can last from a few 
seconds to less than about 1 hour, depending on the number of models in the subgrid. An 
example of - distributions (actually, pj^2 = e~^'^'^T^ with respect to the various stellar pa- 
rameters is presented in Figure [21 



Step 4-' Computation of mean values and uncertainties 



The previous step provides the Xj^ - distributions for each of the six parameters derived through 
the spectroscopic analysis (Teg, log 5, ene, Ct, logQ, and (3). If the absolute visual magnitude 
is provided, the tool automatically determines R, logL, and Mgp (spectroscopic mass) for each 
model, and hence the corresponding Xj. - distributions for these parameters are available as well. 
If also the terminal velocity is provided, a similar computation is performed for the mass-loss 
rate (using Q, R, and Voo)- Finally, one of the IDL modules computes the evolutionary masses 
(Mev) of the models from an interpolation in the (log Tpg, log 5')-plane using the tracks provided 
by a stellar evolution code. Since the - distributions for T^g and log 5 are computed and 
stored in Step 3, the computation of the corresponding distributions for R, logL, Mgp, M^v can 
be easily repeated in a few seconds, in case a different and/or evolutionary tracks need to 
be considered. 
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Figure 2. Example of output provided by the lACOB grid-based tool, for the quantitative 
spectroscopic analysis of the Galactic 09 V star HD 214680. See text for the explanation of 
all the information contained in the figure. Observed spectrum from the lACOB spectroscopic 
database (i? = 46,000). vsmi and Ort were previously determined using a combined Fourier 
transform + goodness of fit technique (see e.g. [14]). In this example, only 5 H/He lines were 
used for the actual analysis. Note that also all other synthetic lines fit the observations perfectly 
at the derived parameters. The analysis was performed in ~ 2 min. See [2] for details on the 
distance, apparent visual magnitude and visual extinction used to determine the absolute visual 
magnitude (Mv). 



These distributions are then used to compute mean values and uncertainties for each parameter 
(taking into account that models above a given threshold in can be discarded). 

The tool also allows to easily investigate possible degeneracies (for example, for stars with 
weak winds the analysis of the optical H/He lines only provides an upper limit for the Q pa- 
rameter, and does not allow to constrain /3), and the contribution of the various parameters to 
the final uncertainty (by fixing one of the free parameters to a certain value and recomputing 
the statistics for the other parameters). 

Step 5: Visualization of the results 



The last step is the creation of a summary plot for better visualization of the results (see 
Figure [2] as an example). The present version of the tool includes: 

• - distributions for the various parameters involved in the analysis (upper-left panels) . 
In those panels, vertical lines indicate the limiting values adopted for the six free parameters 
(i.e., the sub-grid of FASTWIND models for which a is computed), and horizontal lines 
indicate the Px^ value corresponding to (threshold) . The user can select three of the 
models that will be indicated as red, blue, and pink dots (the same colors are also used in 
other parts of the figure). 

• a summary of the input parameters (vsini, Grt, and M^) and the result from the statistics 
for each parameter resulting from the analysis (upper-right part of the figure). 

• a logTcff — logg diagram including evolutionary tracks (e.g., from [11], without rotation, in 
the figure) and the position of all models with X^ < X|, (threshold). 

• a set of panels (lower part of the figure) were the observed and synthetic H/He lines indicated 
by the user are compared. The spectral regions that are used to compute the Xj. for each 
line are indicated in black, while the clipped (or not used) regions are presented in green. 
Note that the plotted lines are not limited to those used in the analysis process (the latter 
are marked with an 'X'). 

2.3. Problems that can be investigated with the lACOB grid-based tool 

The fast performance and versatility of the lACOB grid-based tool not only allows to analyze 
large samples of 0-type stars in a reasonable time, but also to easily investigate the effects of 

• the assumed line-broadening: The common strategy in previous spectroscopic analyses 
of 0-stars was to assume pure rotational line-broadening (in addition to natural, 
thermal, Stark/collisional and microturbulent broadening included in the synthetic spectra). 
However, recent studies have shown that an important extra-broadening contribution 
(commonly quoted as macroturbulent broadening) affects the shape of the line-profiles of this 
type of stars (see e.g. [15], and references therein). How much are the derived parameters 
affected when this extra-broadening is neglected? Can both broadening contributions 
be represented by one 'fake' rotational profile (fjsini), without disturbing the resulting 
parameters? What is the effect of the uncertainty in the derived/assumed broadening? 
Some results from a preliminary investigation of these questions can be found in [16]. 

• the placement of the continuum: The normalization of the spectra of 0-type stars is 
sometimes complicated, especially in the case of echelle spectra. It is commonly argued 
that the derived gravity can be severely affected by the assumed normalization, but, to 
which extent? To investigate this effect one could use the lACOB grid-based tool, applying 
small modifications to the continuum placement. An example of this type of analysis for 
the case of low resolution spectra of B-supergiants can be found in [6]. 



• neglecting/including a variety of different commonly used diagnostic lines: Before computing 
the global (see equation 3), the resulting X| for each line are stored. This way the user 
can easily recompute Xj., discarding some of the initially considered lines. This option 
allows, for example, to investigate the change in the wind-strength parameter when both 
Ha and Hell 4686 lines or only one of them are included in the analysis. 

• clipping part of the lines: Most of the 0-stars are associated to Hii regions. In some 
cases, the stellar spectra can be heavily contaminated by nebular emission lines (mainly 
the cores of the hydrogen lines, but also the Hei lines). This contamination must be 
identified and eliminated from the stellar spectrum to obtain meaningful results from the 
spectroscopic analysis. When the resolution of the spectra is high, the clipped region is 
small in comparison with the total line width; however, even for a moderate resolution an 
important region of the line needs to be eliminated. What is the effect of such a strong 
clipping in the Ha line on the determination of the mass-loss rate? How large is the effect 
of clipping the core of the He i lines in the Tefj determination? 

• the way the statistics is derived from the Xj^ - distributions: As commented in Section 12. 2^ 
this is a critical point in the determination of final best values and uncertainties. 

These are only some examples of things can be investigated with the lACOB grid-based tool. 
Many other tests can be easily performed. In addition, we expect the tool to be of great benefit 
for the analysis of the 0-star samples included in, e.g., the ESO-Gaia, VFTS, OWN, lACOB 
and other similar large surveys. 
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